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Hitting times to spheres of Brownian motions 
with and without drifts 


Yuji Hamana and Hiroyuki Matsumoto 


Abstract. Explicit formulae for the densities of the hrst hitting times 
to the sphere of Brownian motions with drifts are given. We need 
to consider the joint distributions of the hrst hitting times to the 
sphere and the hitting positions of the standard Brownian motion 
and explicit expression for their Laplace transforms are given, which 
are different from the known formulae in the literature and are of 
independnt interest. 
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1. Introduction 

For d ^ 2, we denote hj B = {Bt}t^o a standard d-dimensional Brownian mo¬ 
tion starting from a hxed point x G which is dehned on a probability space 
We throughout assume x 0. Letting v G R'^ be a non-zero con¬ 
stant vector, we consider a Brownian motion B^'^l = with drift v given 

by bI^^ = Bt + vt. It is a very simple fundamental dihusion process, but we 
sometimes encounter difficulty to obtain explicit formulae on it. 

In this paper we consider the hrst hitting times a and of B and B^'"\ 
respectively, to the sphere S!^~^ with radius r > 0 and centered at the origin. 
The main purpose is to give an explict expression for the density of For 
this we need to give an explicit expression for the joint Laplace transform of the 
density of (a, B^^) G (0, oo) x The formula for (a, Rg-) obtained in this article 
is of quite diherent form from the formulae obtained by Aizenman-Simon [1] and 
Wendel [13] . 

The density Pu{t] x), u = being the index, of a has been studied from old 
times. See [Zli and the references therein for the Laplace forms and related top¬ 
ics. Recently, Byczkowski and Ryznar [5], Uchiyama m and the authors of the 
present paper in 0 E] have studied the explicit expressions and the asymptotics 
of the densities themselves and the tail probabilities. 

The density for is expressed in terms of the densities p^{t; xfis (of different 
dimensions). Moreover, using the previous results for ex, we show the asymptotics 
of the tail probabilities for a^'"\ 

Our main results are the following. 
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Theorem 1.1. When d = 2, the density x) for is given hy 

Po\t]x) = 



When d ^3, it is given by 


p^;\t;x) =2"r(z/)e-<"’">-5l"l'* 


oo 

n=0 


Y.(" +")c. 


{v,x) 


\v\ \x\ 


^ i/+n 





a: ” 

n 

Uj>u-\-n 


p^+n{t;x). 


( 1 . 1 ) 


( 1 . 2 ) 


Here {p ^ 0) is the modified Bessel function of the first kind and is the 
Gegenbauer polynomial. 

We refer to Magnus-Oberhettinger-Soni [8] and Watson na abont the special 
fnnctions. We note again that explicit expressions and the asymptotic behavior 
of Puit] x) are known. 

For the asymptotic behavior of the tail probabilities, we show the following. 
To mention the resnlt, we recall (lasin]) that 

Po(t;a:) = ^^^^(l + o(l)) and pfifix) = + o{l)) (i^ > 0) 

holds as t —)■ oo, where L(0) = 2 log and, for z/ > 0, 


m 




Theorem 1.2. Assume |a;| > r. Then, one has 

Pit < < oo) = ^^/o(|n|r)e-<"’"> 

as t ^ oo when d = 2, and 


e-h\Wt 

f(logt)2 


(1 + 0 ( 1 )) 




when d ^ 3. 

This paper is organized as follows. In the next section we give some estimates 
for the modihed Bessel fnnction and the Gegenbaner polynomial. In Section 3 we 
present an explicit form for the Laplace transform of plf \t]x) and, admitting it 
as proved, we give a proof of Theorem 11.11 In Section 4 we show how to compnte 
the Laplace form. In the last Section 5 we prove Theorem 11.21 

We can apply the resnlts in this paper to a stndy on the Wiener sansage of 
the Brownian motion with drift. It will be discnssed in a separate paper. 
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2. Preliminary estimates 


In this section we show some estimates for the modified Bessel function and 
for the Gegenbauer polynomial 
We firstly show an estimate for 


m = 

m=0 


(^/ 2 )M+ 2 m 

r(m + l)r(m + z/ + 1) 


Lemma 2.1. For /r ^ 0 and n ^ 1, one has 

Cn 

Proof. Note that r(p + g) ^ r(p + l)r(g) holds for p ^ 0 and g ^ 1, which can 
be seen from 

r(p + i)r(g) ^ ppip^ q) = p [ x^~^dx = i, p > o. 

r(p + g) Jo 

Then we have 

r-Mj ('£') < Jl_ ^ _ 

? (r(m + l))2r(p + n + l) 

^ c i^/^r y 

~ 2f^+'^T{p + n + ml ) ’ 

which shows fl2.ip . □ 


Next we give an estimate for the Gegenbauer polynomial G^. When z/ > 0, if 
is given by 

1 FI m I 


TIf) 

' ' m=0 


m\{n — 2m)! 


which is characterized by the relation 

CXD 

(i-2(f+0“'' = I]cKr. 

72=0 

When F = 0, Go(0 = 1 and, when rz, ^ 1, G° is given by 

r(n — m) 


[n/2] 


m=0 


r(m + l)r(n — 2m + 1) 


(20 


72 — 2772 


Lemma 2.2. For a G R with |q;| ^ 1, f d 0 and n'Pl, one has 

4”r(F + n) 




n\ 


( 2 . 2 ) 


where po = 1 o^'^d py = (r(F)) ^ for f > 0. 
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Proof. When i/ = 0, since (^^(cos^) = - cos(n0), we have 


n n\ 


When z/ > 0, we have 


< 


[n/2] 

E 


1 2"' ^™'r(z/+ n — m) 


r(z/ + n) r(z/) r(z/ + n)m!(n — 2m)! 


If 1 ^ m ^ [n/2], it holds that 


r(z/ + n — m) 


< 


r(z/ + n)(n — 2m)! (n — l)(n — 2) ■ • ■ (n — m) • (n — 2m)! 


(n — m — 1) • ■ • (n — (2m — 1)) ^ n 


(n — 1)! 


n! 


Hence we get 


r(z/ + n) “ r(z/) ^ m! n! 

^ 2” ^ 1 /n\”^ ^2" t < 

“ r(z/)n! ni! v4/ ~ r(n)n!^ “ r(z/)n! 


because e4 < 2. □ 


3. Laplace transforms and proof of Theorem 11.11 

We hrst reduce the computation for to that for the joint distribution of 
(a, Ba-), the hrst hitting time to the sphere and the hitting position of the standard 
Brownian motion. 

By the Cameron-Martin theorem, we easily see 




and 


^[g-AaW] _ 


where B denotes the expectation with respect to P. Moreover, letting = 
o'{Bs, s ^ t}, we see by the strong Markov property of Brownian motion 


and 




/•CXJ 

^[g-AFD] _ yg-{^,x> / 


From this identity we obtain the following 
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Proposition 3.1. For any A > 0, one has 

^[g-AaO)] ^ 

Proof. We have shown 


/ e-^^dt / = s]p^{s; x)ds, 

Jo Jo 

where Pu{s] x) is the density of a (see Sect. 1). Changing the order of integrations, 
we obtain 

poo poo 

^[g-A<xO)] _ / E[e<^’^‘^>-^l’'l''^|cr = sK(s;x)ds / e-^^dt 


= g-h.a:}^|'g{i),Ra>-(A+2hr)o-j_ q 


For the right hand side of fl3.ll) . we show the following explicit expression. 
Denoting by the modihed Bessel fnnction of the second kind (the Macdonald 
fnnction), we dehne the fnnction (p ^ 0) by 



a;Kv/2a + ) 

A'p()(\/2ATyp) 


if .^ > 7] > 0 


and 


vO.A 


i^,v) 


I,i^^2X + \v\^) 

4(r/^2A + |n|2) 


if p > ^ > 0. 


Since is decreasing and is increasing on (0, cxd), ^ 1. 

Proposition 3.2. Let A > 0. When d = 2, one has 


n=l 


where a = rr^. When d > 3, one has 
hl-hl — ’ 


CXD 

= 2^r(z/) y (n + n)Ci;(a) 


n=0 


Iu+n{\v\r 
(|n| • |a:|) 


_ vfofA 


(|a;|,r) 


Combining this proposition with fl3.ip . we obtain 
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when d = 2, and 



n=0 


Iu+ni\v\r) 
(|n| • \x\Y 


^v-\-n 


{\x\,r). 


when d ^ 3. 

We postpone a proof of Propositon 13.21 to the next section and give a proof of 
Theorem o It is well known (cf. [7]) that, the density p^{s;x) of the hrst 
hitting time of a Bessel process with index p starting from |a;|, is characterized 
by 

E[e-^^Y= j e"^*PM(s;x)ds = —Z®’^(|a:|,r). (3.2) 

Jo 

We use the same notation for the density since our main concern is on the special 
case where the index /r is a half integer, and there is no fear of confusion. 

To prove Theorem 11.11 we compute the Laplace transform of the right hand 
side of fll.ip . fll. 2 p by changing the order of the integrations and the inhnite sums. 
Using the estimates given in the previous section, we have for z/ ^ 0 


CXD 

E 

n=l 


(z/ + n)|C';((a)|/^+„(|n|r 

OO 

E 



X 


VrpV-\-n 




< JYlLpMr Y + n)T{u + n){2\v\r)- 


(3.3) 


2^ X 


n=l 


n\V{v + n + 1) 


^ u-\-n 


(|x|,r). 


Since E 


1 , the above is bounded by 


PyT’' 

2 ^ 1 x 1 


Mr 


E 

n=0 


( 2 |x|r) 


n\ 


PM 

2 ^ 1 x 1 




Hence, we may apply Fubini’s theorem and see, from fl3.2p . that the Laplace 
transforms of the right hand sides of fll.ll) and fll. 2 p are equal to those of x) 
in both cases. 

We have now shown Theorem 11.11 admitting Proposition 13.21 as proved. 


4. Proof of Proposition 13.21 


In order to prove Proposition l3.21 we use the skew-product representation of Brow¬ 
nian motions. Let R = be a d-dimensional Bessel process (with index 

u = 4^) cind 9 = {9t\tzo be a Brownian motion on the unit sphere S'^~^ = Sf~^ 
with 6^0 = ifj, and assume that R and 9 are independent. Recall that, embedding 

gd-i 

we can realize 0 as a solution of a stochastic differential equation, 
which is so-called Stroock’s representation of a spherical Brownian motion. 

Set St = jl{Rs)~‘^ds. Then, {Rt9st }t^o is a d-dimensional Brownian motion. 
Hence, we have 




'OO 


G df, F, G dw), 

|a:| 
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where denotes the expectation with respect to the probability law of 0 starting 
from 6 * 0 , -Pi/ p| is the probability law of {Rt} and r is the hrst hitting time to r of 

It is known (cf. [2] p.407) that 


Ixl '^K 




( I ^ I ) 


r-i^K 






if \x\ > r 


and 


|x| 




y'i/2|^2(l^l'^2tt) 


r-^J 




if |x| < r, 


where -Pi/,|a;| is the expectatation with respect to Pi,,\x\- 

We obtain Proposition 13.21 if we show the following. We can justify the change 
of order of the integration and the inhnite sum by the same way as fl3.3p . 
Proposition 4.1. Let ^ > 0. Then, when d = 2, one has 


= 7„(|„|r) + Y,nCl(a)e-^’'nMv\() (4.1) 

71=1 

and, when 3, 

= 2’'T(i') ( 4 . 2 ) 


Kl 


n=0 


We see from this proposition that the Gegenbauer polynomial comes into our 
story through the following formula (cf. [H p.227]): for a G R, ^ > 0,/r > 0, 


e“S = 2^(1,) + n)C;(a)r'‘/„+n(0. 


(4.3) 


n=0 


We hrst show that fu{t,0 = s^-fishes 

df^ d-l^dfy , 1 


= + t>o, e>o, 


(4.4) 


(4.5) 


dt 2^ 2 

together with the boundary conditions 

/i/(0,O = fu{t,Q) = 1, ^(^,0) = (T,0o)e"'^^ 

For this purpose, we recall Stroock’s representation of sperical Brownian mo¬ 
tion (cf. [in])- 0 may be realized as a solution of the stochastic differential 
equation based on a d-dimensional Brownian motion {tCs = {wl,Wg,...,w^)}s^o 
which is given by 

d 

i=i 


i = 1,2, ...,d. 

















8 


Y.HAMANA AND H.MATSUMOTO 


Then, from a strihgtforward computation by using Ito’s formula, we can show 
(14.411 ■ It is easy to see (14.511 . 

For simplicity we set /? = d — 1 and consider the function g^, given by 


9At,0 = T^)- 


Then is a smooth function which satishes 

dgt. ^2d^9u , ^2 n n 

dt ^ de ^ ^ ^ > 0 ’ ^ > 0 ’ 


(4.6) 


and 


9u{0,O = e'"^, 9u{t,0) = l, 


d9v 


(t, 0) = ae 




(4.7) 


If u{t,^) = e ^*0(0 satishes (14.6p . we should have 


eV"(o+/9e0'(o - + mo = o. 


The system of the fundamental solutions of this second order differential equation 
is given by and where u = ^ = ^. For the 

function 0 to be smooth at ^ = 0, we should choose Moreover, 

n = \/A + z/2 — u should be a non-negative integer and A = n{n -|- 2iy). 

The following lemma is easily shown and we omit the proof. 

Lemma 4.2. (1) The function (pu,n{.0 = ^~’^O+n{.0 satisfies 

om(o+mu,nio - OtuAo = '^ > o- 

(2) One has 

= 2<^+^T{v + 2) ^ ° 

The following proposition immediately implies Proposition 14.11 
Proposition 4.3. When d = 2, one has 

OO 

9 o(t,i) = /„(0 + C“(a)e-"’‘/„(J), t ^ 0, ( ^ 0 (4.8) 

n=l 

and, when d ^ 3, 

OO 

gAt, 0 = 2U(0 + n)C;;(a)e-”‘"+="y-V„+„(0. (4.9) 

n=0 

Proof. First of all we note that the sum on the right hand sides of (14.8p and 
(14.911 are absolutely convergent at each {t,f), which is seen from (12.ip and (12.2p 
in a similar way to (13.3p . 
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Letting be the function defined in Lemma [4.21 we set 

OD 

K{t\ 0 = !]('' + n)C,';(a)e-”l”+"-V,n( 0 - 

n=l 

We have shown that the sum on the right hand side is absolutely convergent. 
Moreover, noting 

Tl 

we see, in a similar way to fl3.3p . that 


El" + and + n)G(n)e-"'”+"">V'',„(0 


n=l 


n=l 


converge uniformly on compact sets in ^ G (0, cx)) and are equal to and 

respectively. 

Next we look at hu{t,^) as a function in f > 0. By fl2.ip and fl2.2p we have 

CXD 

E(<^ + n)\C:(a)\n(n + 2n)e-”'“+=”V,n«) 

n=l 

/ A^Tiv + n) , , r 

- ^ ' n\ ^ ^2^+"r i/ + n + l 

n=l ^ ' 

^ ^ . f V , V (2n + 1)(20' 

2- (n-1)! ^ 


n=l 


n=l 


{n — 1)! 


= , P^(2l/+l)e 3£ 

2U-2 2^-1 


and we may differetiate term by term to obtain 
d 

hu[t,^) = - 

n=l 


dt 


hu{t,0 = + n)Cni(^)e + 2z/)v?^,„(0. 


Combining the identities above, we see that the function gu{t,^) given by fl4.8p 
and 04.91) satishes 04.61) . 

The boundary condition 04.7p in the case of d ^ 3 may be checked by 04.3p 
and the fact, C'o(a) = 1 and Ci{a) = 2z/a (cf. [H p.218]). 

For a check when d = 2, we rewrite 04. 3 p as 


e”« = 2'‘r(^ + l)r‘‘IA() + + n)CS;{a)r>‘I,+n{i) 

n=l 


= 2'‘r(f, + i){r%K) + f^ip + 

'n — 1 ^ 
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which holds for any fi > 0. Note —)■ C°(a) as /i 0. Then, by nsing fl2.ip 

and fl2.2p . we can show that we may apply the dominated convergence theorem 
and obtain 

OO 

n=l 

which is exactly 5'o(0,0 = 

Another bonndary condition gQ{t, 0) = 1 follows from /o(0) = 1 and the other 
one 0) = ae“* does from the formnla Ci{a) = a. 

We have now completed the proof of Proposition 14.31 
Remark 4.4. The fnnction z~'^ may be regarded as a holomorphic fnnc- 
tion on C. Hence, by nsing (EHi cLiici fl2.2p , WG CcLii show tlicit tho furictioiis 

OO 

n=l 

are holomorphic in z G C. From this we obtain the Fonrier-Laplace transform of 
the joint distribntion of (a, Ro-). For example, when d ^ 3, we can show 


= 2"F(0 ^ + n)Ci;(a) 


n=0 


Ju+n{\v\r) 
(|n| • \x\Y 


7ig),\ 

-^u-\-rL 


(|a;|,r). 


where is the usnal Bessel fnnction. 

5. Proof of Theorem 11.21 


We set 

/ OO 

x)ds. 

In order to apply Theorem 11.11 we need to change the order of the integration 
and the snmmation. 

For this pnrpose, we note from fl2.ll) and fl2.2l) 


+ n)\C^Yo)\- 


u-\-n 


(|n|r)|a;| 


n=l 

OO 

n=l 




e 2 1^1 ®pj,+„(s;a;)ds 


4”F(z/ + n) (|n|r)"'el’'b |a;|"' _iu,i 2 

n\ 2^+"F(z/+ n + 1) r” 


^\v\H 


~ 2’^ ^ n\ ~ 2^^ 

n=l 

Then we can apply Fnbini’s theorem and obtain from Theorem 11.11 

/ OO 

e~^^''^^"po{s;x)ds + fo{t) 
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when (i = 2, and when d ^ 3 


+ 1 ) 


Iu{\v\r) 


'X) 




where the second terms on the right hand sides are given by 
Uit) =e-<"’">2"r(n+l) 

\x\^Iy+n{\v\r) 


X 


+ n)C'^ 




n=l 


\y\Vrj^V + n 


e 2 l"'l *p^+„(s;a:)ds. 


At hrst we prove the theorem when z/ > 0. We have 

Pv{t]x) = ^^(l + o(l)) 

and, by L’Hospital’s rule, 

x)ds = + o(l)). (5.1) 

Note that this identity holds for any n > 0. 

In order to show that fu{t) is negligible when z/ > 0, we use the argument in 
Sect.2 of [B]. Then we obtain 

^ < oo) 

where E^^n,\x\ denotes the expectation with respect to the probability law of the 
Bessel process {Rt}t^o with index u + n and starting from |a:|. 

Using the explicit expression of the transition density of the Bessel process, 
we obtain 






< 


(2f)'^+" 

1 

i^+n 


(2i) 


_j^ 

e 2 t 


nr 
e 2 t 


X 


E^r{v+ n + m + l)(2ty 




\X\ 


^■^0 r(i^ + n + l)r(m + l)( 2 f)' 


r(z/ + n + l)(2t)'^+”' 

Hence, by using fl2.ip and 02.21) again, we obtain for n ^ 1 and f ^ 1 


4.V+1 


el'"f‘(v + n)\C:(a)\ 


X 

"4+n(|n|r) 


n 

UrpU-\-n 


< p-^M, (ri-ri)“ 

“ 4^ n!r(z/ + n + 1) 


e-|H= 


'p^+n{s]x)ds 
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The quantity on the right hand side is independent of t ^ 1 and is summable in 
n. Therefore, since fld.ip holds when we replace z/ by z/ + n, we can apply the 
dominated convergence theorem and obtain 

lim = 0. 

t^OO 


When d = 2, we have 

+ L(0) = 21og^, 

and 

x)ds = + o(l)). 

By the same way as in the case of z/ > 0, we can show by the dominated conver¬ 
gence theorem 

lim f(logt)2e^l"l'VoW = 0 

t—^OO 

and we obtain the assertion of Theorem 11.21 also in this case. 

Remark 5.1. The estimate for the tail probability P{t < a < oo) has been 
hrstly given in Byczkowski and Ryznar [3]. We need an explicit upper bound 
here. 
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